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MOTION OF A CYLINDER UNDER THE SURFACE OF A HEAVY FLUID" 


By L. N. Sretensky 


1. INTRODUCTION 


The present work on the theory of the motion of a solid body in a 
fluid having a free surface consists of two parts. 


In the first part (sections 2 to 8), general equations are given 
for the determination of the flow of a heavy fluid of infinite depth 
about a submerged circular cylinder. The problem of the motion of a 
cylinder under the surface of a heavy fluid presents considerable dif- 
ficulties in its solution. These difficulties were first pointed out 
by Kelvin. A solution is given herein for the simplest part of the 
problem of Kelvin, namely, setting up the equations of the problem and 
obtaining certain approximations of their solution. The approximate 
solution obtained replaces the moving circular cylinder by a certain 
vortex. 


T. H. Havelock (reference 1) in a recent paper considers the prob- 
lem of Kelvin under the same general assumptions as are herein con- 
sidered, but gives a more advanced approximate solution. 


In the second part of the paper (sections 9 and 10); the cylinder 
is replaced by a dipole of a certain strength, and an equation is set 
up for the computation of the wave resistance of a circular cylinder 
moving in a fluid of finite depth. 


2. DERIVATION OF BOUNDARY CONDITIONS OF PROBLEM 


In order to study the problem of the motion of a cylinder under the 
surface of a fluid, a system of Cartesian coordinates XOY is intro- 
duced. The OX-axis is placed along the undisturbed surface of the 
liquid with its positive direction coinciding with the velocity c of 


"uDvizhenie tsilindra pod poverkhnostyu tyazheloi zhidkosti." 
CAHI Rep. No. 346, 1938, pp. 3-27. 
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the flow at infinity, and the OY-axis is taken vertically upward. The 
radius of the cylinder is denoted by a; and the ordinate of its center 
by h. The origin of the XOY system is chosen directly above the 
center of the cross section of the cylinder. 


The fluid moves with the velocity potential & and the stream 
function v. The disturbance (due to the presence of the cylinder) of 
the horizontal projection of the fundamental flow velocity c is 
denoted by u; the disturbance of the vertical projection, by v. If 
c is the horizontal velocity of the flow at infinity, the following 
expression can be obtained for the projection of the total velocity v 
of the particles of the fluid: 


u + c,v (1) 


The eguation of Bernoulli can now be written, if the motion is assumed 
to be steady: 


P Li gi 
ania Sek N (2) 


where O is the mass density of the fluid. 


Along the free surface the pressure has the constant value Po- 
Eguation (2) is applied to the particles of the free surface lying far 


ahead of the cylinder, that is, to those for which x= - æ. For 
xs: - e, V < c and y= 0; therefore, 

Po 

o. -l še 

5 = C 5 c 


Equation (2) can then be written 


P-Po P ase 32 
Day - L (v-e?) 
But 
ve = (urc)? + vč 
hence 
P-P 
-= = - gy - > [v2 + v? + 2cu] (3) 


The increments u and v are now assumed to be so small that their 
squares may be neglected. With this assumption, the accurate equa- 
tion (3) becpmes the approximate equation 


O ei (rusiji 
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When this equation is applied to the free surface, 


ki aw ake: AN 


This equation can be used to determine the shape of the disturbed sur- 
face of the fluid. 


The fact that the velocity of the particles moving along the surface 
is directed along the tangent to the surface is now taken into 
consideration: 


V 


dy _ 
ax ute 


This relation may be represented in the form 


a Z A 
dx c Ge TAS 


If the terms of the second and higher order smallness are rejected, 


dy v 
ax" c (5) 


Eliminating the ordinate y by combining equations (4) and (5) yields 


du 


= Kalla = v (6) 


By introducing the potential (x,y) and the stream function ¥(x,y) 
of the disturbed velocities, this condition may be given another form: 


de _ Ov 

ee Ox 3 

v=- (7) 
dy ox 


The functions © and VY are harmonic conjugate functions. Substi- 
tuting expressions (7) in condition (6) yields 


di = & Oy (8) 
Ox Oy ci Ox 


A second simplification of the problem is now introduced by requiring 
that y be replaced by zero. Eguation (8) will then be satisfied 





4 NACA TM 1335 


along the axis of abscissas. Integrating along this axis gives 
dv g 
= —V (9) 
dy 2 


where the constant of integration is zero from the consideration that 
the zero value of the stream function Y was ascribed to the surface, 
and moreover, for x = -, u= - dW/dy = 0. 

Condition (9) is the first boundary condition of the problem. The 
conditions on the surface of the cylinder will now be discussed. Let 


the cylinder be washed by the streamline y =a. Since the relation 
between Y and v is given by 


V - -cy t v 
the condition Y= a on the cylinder will take the form 
V= a+ cy 


where a> 0. 


The determination of the motion of the fluid thus depends on the 
integration of the Laplace equation 


ov 
ay = i v for y = O (1) 
V < a t cy for x? + (y+h)? = a? (II) 


5. TRANSFORMATION OF FLOW REGION INTO CIRCULAR RING 


In order to obtain the integral of the Laplace equation correspond- 
ing to conditions (1) and (IT), a conformal mapping of the complex var- 
iable plane z = x + iz on the plane € = &+ in is carried out by 
setting 





Kan EET 


Sa 


er ox 
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ts z (1) 
where 


A= - Sh? - 22 <0 (2) 


When the complex variable z is varied along the OX-axis from -œ to 
te the point € describes, starting from the point € = 1, the circle 


Eö + ne = 1 in the clockwise direction. 


If z=- hi + ae} is set next in formula (1) and the angle y 
is varied from zero to 2x, the surface of the cylinder will correspond 


to the circle ze + nê = p? of the plane of the complex variable €. 
The radius p of this circle is determined in terms of a and h by 
the formula 


_ he /hi-af 
P a (3) 
or 
htA 
=. 3! 
p = (34) 


The number p is evidently less than one; h > a. 


The region occupied by the flow is thus transformed into the space 
enclosed between the two concentric circles ig] = 1 and | EI = p. 

The numbers A and p introduced previously have simple geometri- 
cal meanings: A is the length of the tangent drawn from the origin of 
coordinates to the cylinder, and p is the tangent of the quarter 
angle & subtended by the cylinder at the origin. 


From formula (1) z is determined: 


z = Ad Tt (4) 


and from this value are found the values of x and y in terms of the 
argument w and the modulus x of the complex number €: 
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2Xr sin W 
ss n 
1 - 2r cos O +r 


x a (1-rê) (5) 


l - 2r cos O + rê 


r sin w= - << 
x?" + (y+X) (51) 
ROM 
r cos O = TOIMIN 
x" + (y+) 


The transformation of boundary conditions (I) and (II) of the problem 
follow. 


Condition (II) will be considered first. Formula (5) shows that 
this condition may be written as 


2 
[x] osat a e (TI') 
zob l - 29 cos o + p 
In transforming condition (I) it is noted that from formula (5) 
there follows for' y = O 
Adr 
1 cos @ 
from which condition (I) assumes the following form: 
Bi et] (rt) 
or l - cos W rel 


where 


Conditions (II') and (I') may be given another form by introducing the 
complex stream function: 


w= + iv (6) 
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The conditions then assume the following forms: 


mag. [ziv + (1-6)? zl o for |= (111) 
Imag.(w-cz) =a for |b] =p (TI'') 


The function w is holomorphic at all points of the ring p < |t]< 1, 
but it will be assumed that the flow around the cylinder has a definite 
circulation. The function w will therefore not be single-valued in 
the ring under consideration. 


In order to obtain the reguired function w(t), the form of the fol- 
lowing infinite series will be used: 


w=qin€+ D pit (7) 


with the undetermined coefficients Pm and q. 


4. COMPUTATION OF COEFFICIENTS OF SERIES 
REPRESENTING FUNCTION w(¢) 


First considered are the relations obtained when condition (II'') 
is required to be satisfied by the function w(¢). 


In the notation 


Pm = Am + Bni 
(8) 


the coefficients An and Bn are real numbers, and T is the circu- 
lation of the flow about the cylinder. 


Condition (II'') then easily transforms into 
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+o 


5 2 
+ ln p + > Pu, cos mp + A sin mp) - McR(1+26426 + ...) =a 


from which the following relations are obtained: 


ŽLinp tB < k ta 


2x 

ot + — = 2hep™ 
p mel, 2,3... 
-m 

EN S 
p 


Condition (I'') is now considered. When the following formulas are 
used 


Imag. [ee] = ZuBo + 2u 5 (Bn*B.m) COS mo + (An-A m) sin mp 


m=1 
tag. [(1-0)7 s] =- L + Imag. = € + +) + 
+0 
Imag. J (an - 2mpy, + (1-21) (N 
11=- co 


no difficulty is encountered in representing the condition 
Tuag. | Zu + (1-6 )E s] eni 


in the form of the following trigonometric series: 


00 


2uBo + 2u 2 (atan) cos m + (A,-A p) sin w) = 


+o 


+ E cos % + Je (edn - 2mB, t (1). cos mb + 


m=-s 


AIF 


(rana - 2mA, + (2-1) 8 sin æ) =0 


This equation gives a series of relations between the coefficients A, 


and Bm: 
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By = Bo = T - 2uBo (11) 
W(B,+B 3) = (B,-B 1) - (B2-B.2) - Z (12) 
H(BntB m) = 1(Bn-B.m) - ta (Bn+1-B-(m+1)) g 
TL (Ba-1-B-(m-1)) (m = 2, 4, 6. . .) 
| (13) 
| 14 (An-A m) = (Ath) - BEE (Ansa tA (ml) - 
| == (Am-1+A (m-1)) (m = 1,2,3,4,5. . .) 
(14) 


The last of these recurrent relations will be discussed first. From 
recurrent relation (10) 


am 
Alm =P Am 


as a result of which relation (14) may be written 


m(1+p za An + 


ma GRE NN] ja [14(1-02) 


2m+2 


mtl 


5 O (m = 1,2,3. . .) 


(1+ọ 


)Am+1 
For briefness, 


> (1+90) A, = km 


a E ne 
m 2m m 


1+p 
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With these notations the obtained relation can be rewritten in the fol- 
lowing form: 


Kn-1 + Smim + Km-i = 9 


(130834 x 4) 


(16) 


The coefficient Ao is a certain finite number so that ko = 0. 


Relation (13) is now considered. Eliminating from this relation 
the numbers B with negative indices and making use of relation (10) 
yielä 


E (uoaa, 4401-02) -alaton Jan + 


m+1 2m+2 _ 2mim-l , m+l 2 
= (Me Ce ae SES 


For brevity, 


z (10) By = x, (17) 


is written to obtain for Xxx the following relation: 


m-1 m+1 
Xm-1 t Smm + "adi 7 a ar p? -m bi (18) 
(m = 2, 3, 4, 5, ) 
If 
=. 21 
Xo = p (19) 


relation (18) for m = 1 gives condition (12). 


The system of recurrent relations (16) and (18) is thus obtained 
for the coefficients 


Ag Big} Bg ae is 
Bi, Bo, Bz, . . . 


To these relations there must be added relations (9) ana (11). 
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D 
4 


5. INVESTIGATION OF RECURRENT RELATIONS BETWEEN 
i NUMBERS ky 


The relation 


== 


Km-1 + Smkm + Km] = 9 





























(16) 
ja (m = 1; 2, 3, . œ ;) 
i may be used to set up two functional relations between the two 
v functions 
oo 
F(¢) = Ant” 
m= 
` (20) 
o(£) = >t 
m=1 
of the complex variable €. 
First of all 
k. m 
P(E} = JEST S 
m=1 1+p°™ k 
or 
F(t) = | = Gust sat 
Om=1 1+p°™ 
The new function S(t) is introduced 
km m 
S(t) = ` zn 5 (21) 
mel 1+p 
from which 
6 
F(t)=2| s(t)-% (22) 
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The function S(t) is connected with the function o(t) by the fol- 
lowing obvious relation: 


S(t) + s(p?t) = o(t) (23) 


Together with this relation there is still another obtained from the 
following considerations. Replacing Kn on the right side of equa- 
tion (20) by 


= Sm-išm-i ~ Km-2 
yields 
o(f) = - > Sm-1km- 16" - D km 26 + kil 
m=? = 
or 
(1+6)o(t) = ke - € 2 One E as (24) 


Transforming the infinite sum on the right side gives 


s je 2m-2 m-l 
2 m-1Xm-19 l= 2p 2 i Emin 7 25(£) 
OT J t s 
N __ 2m-2 
3 En-ikn-i6 = ču ; 2 mani Ene ši a 80 


The function under the integral may be expressed as follows in terms 
of the functions o(t) and S(t): 


z 2m- 
a = ft! = 28(t) - o(t) 


me? l+p2m-2 
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hence 


sm-lkm iti = 2 I 25(€)-0(6) ae _ 
J 1km-16 H 5 E € - 20(t) 


Substituting in eauation (24) the obtained value of the infinite sum 
produces a new relation between the functions S(t) and o(€): 


6 
-t)7oft) = s 28(6)-o(£) 
(1-0)50(€) = iyo - aug [ŠE at 


Rewriting this relation in the differential form, 
a |(1-6)"o(6)]_ _ 9, 28(£)-0(t) 
a [anis]. ar < (a 


and adding to this differential equation the functional equation 


S(t) + Sloĉt) = a(t) (23) 


furnish the two necessary equations for the determination of the func- 
tions o(t) and S(t). The function S(t) being determined from these 
equations, the function F(() may be obtained with the aid of rela- 
tion (22) or the relation 


ut 


2 
rit) = fm - SEE oto] + oft) at (22") 


2u o 


The function 


y(t) = 2) ae 


may be represented through the function F(t) as follows: 
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n (6) = #() 


The complex stream function (7) is represented thus: 


W = w + iw (26) 
Setting NE 
me 2 Ab" (27) 
m=- © 
and 
Wo = £ In € + 3 Bao (28) 
m=-0 


the function w,(t) may be constructed from the function F(t) by the 
formula 


v (€) = Ay + F(t) + r (e) (27') 


or 
2m 
m 
aED a(t i 
m=1 c 
In the variable z this function may be written 


m(x) = Ay + J, An « (2-14) Mp a eri) 
m=1 


(224, 2)0 


The function Wo must now be investigated. 
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6. INVESTIGATION OF RECURRENT RELATIONS BETWEEN NUMBERS xm 


In the previously obtained recurrent relation 


m-l  mtl 
Uni + Smm + %m41 = pepali, =o pe - m - J 


(29) 
% = (m=1, 2, 3,.. .) 





the three functions of the complex variable E may be considered: 


Qo 


Ba 
N 
B 





F'(€) = 
m=1 

o'(£) = x ai? 
m=1 

S'(t) = = s 
mel lipo™ 


Between these functions there exist the following two easily obtained 


relations: 
S'(¢) + s'(p?t) = a(t) (30) 
E 
' = 2 1 ja ag 
F'(0)=2] sv(t) (31) 


Equation (30) gives the first relation between the functions S'(() 
and o'(t); for obtaining the second relation, the recurrent relations 
(29) are used as follows: 


(6) = x Sm-1%m-16 - Ja *m-26" tuje + L(G) 
m=2 m=2 
OT = 
L(€) = 2%c J = + Eu o -m- desa (32) 
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Transforming this relation gives 


(1+t7)o'(t) = - € 5 aa a = xot? +x + L(t) 
m= 


s t 
a, saata ai = | AT) ag - 201 (8) 


m=2 


(by the computations of the preceding section). Hence 


c 
(1-t)20'(t) = - os | PEE) CG) ge. xot? + xat + LE) 
from which is obtained 
2 


The function L(t) may be represented in finite form: 


1 (6) _ pi(t-1)8 p pečt 
0 2(1-p%t)@ 2 1-p%t 





whence 
4 ar), prp (ET pi 
8 (1-86) (1-0)? 
a L(č) 


Substituting this value of the derivative < ica in eguation (33) 
yields ag 


2 
d 1- k t g! 
2 [1200000] =. ax SG. ry 


(34) 





(1-p#¢)5 (-078)4 


ep? fae Ja eee | 
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This equation together with equation (30) may serve for the determina- 
tion of the functions o'(€) and S'(t). The function S'(t) being 
found, F'(¢) may be computed by formula (31). 


The function wo of the preceding section may be written as follows: 


so 


Wo(t) = By + Ze: in € + x Batt + È B.,,6"" (35) 


mel 











The second sum may be transformed into 


J Bmt ™ = - 9 2) s J pagn 
m=1 





m=1 m=1 
or 
: -m 2 2cr 2 
B- = - FUE ) p 
2 = G t-p? 
Thus 
E T 2rcp2 t ) 1 p? 
a a Ks NS - F E (36) 
Or 
2 = Va 
welt) = % +E ing + Bees) agro) (37) 
C-p m=1 6 


In the variable z this expression may be written 


z-li , 2rcp“(z+%) 


+ 
ZtAL (1-p2)z-Xi(149%) 





poe ae 
we(z) = Bo + z m 


(38) 


Yy Ba- (Ga) Ds a " 
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7. FORMULA FOR COMPLEX STREAM FUNCTION 


On the basis of formulas (27') and (36), the expression for the 
function w(t) may now be stated 


w(€) = (AptiB) + Hin js + Bor i + 


(39) 
F(t) se) F(t) s (B-) 
+ + - = 
c va L c E 
Or 
00 oo m 
f ri 2rcp@ — fe? 
w(t) = (AptiBo) + 5> In € + = i + X Pace + x zdi 
RA b-p? m=1 m=1 E 
(39') 
The number Pa is the conjugate complex of the number 
Pm = Am + iB, 
Equations (9) and (11) permit the determination of the coef- 
ficient Bo and the number a: 
T 
By zač n 
Ho lea 
(40) 
Bak 
ae US EEA 
a= On inp + i-2u Ac 


The number A, remains arbitrary, and the circulation I likewise 
remains arbitrary. 


Formula (39') may be given a more symmetrical form by introducing 
in place of the coefficient p, new coefficients Tn and Yi by the 
formulas 


rn = Pm ~ tle 
(41) 


Pm + ide 


5 
HI 
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to obtain 
; Tä acp? , Xet 
w(t) = (45 +189) + ze In € + E + it + 
00 œ 2 m (39") 
2 rat +) Fall) 
m=1 m=1 


8. APPROXIMATE SOLUTION OF PROBLEM OF DEEPLY SUBMERGED CYLINDER 


The problem of determining the flow of a heavy fluid about a cir- 
cular cylinder is thus reduced to the solution of the two systems of 
the functional equations (23), (25), (30), and (34). 


Since it is not possible to solve these equations without the aid 
of infinite series or in a finite combination of elementary functions, 
an approximate treatment of the equations of the problem is proposed. 
For this purpose the number p? must be considered. In the following 
table are given various values of the ratio a/h and the corresponding 
values of p and pt. 


609|0.0718 |0.005151 





This table shows that for a ratio a/h less than 0.342 # 1/3, the 
number pt does not exceed 0.001, which justifies the rejection of all 
powers of p? starting with the second in considering the motion of 
the stream for a/h < 1/3. The results here obtained are in somewhat 
complicated form; therefore only formulas which are suitable for the 
condition at which it is permissible to reject the components with p? 
are presented. The preceding table shows that this can be done by 
starting, for example, with a/h < 1/3. 


By rejecting the terms with p? equations (23) and (25) may be 
rewritten as follows: 


S(t) = a(t) 


a 
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from which 





The function F(t) is defined 


from which 


Equations (30) and (34) are now considered. These, with the preceding 
approximation, may be written as 


S'(t) = o'(t) 


a jem ot] «oy 1) E 


at an 


from which 
2 S aut 
ot(€) = ETÄ 1-6 x] - = 5 ae dt 


The function F'(t) is now obtained: 


6 
e) galeb) LT 1)(1-6)"o'(€) 
roe] eee eS [Gita J 
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The complex stream function may now be found. For this purpose for- 
mula (39) is used to obtain 


= 
k m k, hei -1 
1 Ti 1 1 

zz P= oe a ei. + 
w(t) = (AotiBo) pe È In E + 7 TA t n e 


(41) 


ee 


The constant By) may be determined by formula (40). Since B = 2x 
(by relation (17)), 


O` 1-2u 


In order to determine the remaining constants kj and xj the con- 


dition of the absence of disturbance of the flow ahead of the cylinder 
is used. The original variable z in formula (41) is used: 


ba. Z= Xi z-ki 
(£) = (45 + 2 TTT 120) - Ku” ZEAL = ZE na" 





ky jiä u c r c? 


dz 
n e se oe e nč (E VITA 


Integrating by parts, 





gz, 
- Si 
k i ; i  kjHuji Z 
N = , lih, DE jn eM LE He < 
w(t) = (4 pa + Boi - ci z Ta n 7 eHe + 
- &Z i s ai 
2 c 
Ti oC Soo Ge 
n Xi zhi 
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E 1 ot Ji z 
es 
— dz + 
ny IN M = 
œ z EZ . 
E cos sš + 2A sin g 2gh 2 = 
2 ae c2 ec 
= g Dar dč + xie + SY az 
Hence 
a Ti Ti z-Xi 
v(6) = (80 + + Boi - n), pe lat 
2gA co N 
REDI S : Ecos SŠ + 2A sin & - Bia 
a eh. Te" + li = déle < + 
ši O 4222 
n poi 
Ti ce ef 
— e dz 
x ZEAL 


The condition at which the function Y becomes zero for y < 0 and 
x -e will be satisfied if the brackets are equated to zero: 


ky — oe uje 


a E cos = + 2X sin si 
ur -u c c 

a a aa aN 

= b ANGLE 


and in addition the imaginary part of the parentheses is equated to zero: 


svile 
nu 


The number Ap is arbitrary; the real part of the terms in the paren- 
theses is also equated to zero to give 


In this way all the constants have been determined and the complex 
stream function of the total flow may be written 
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-Zi (" Z l 
_ r Z-hi T c e 
Poa Simo aLt  <[ iza (42) 


Hence for the case here considered of the deep submersion of the cylinder, 
the flow consists principally of the circulation of the velocity about 
the cylinder. The sources of 1 n strengths entering the general 
formula (39") located at the point z= Xi begin to be effective only 

at small depths of submersion of the éylinder. 


The vortex representing the cylinder is situated at the point 
z= Xi, this point being somewhat displaced with reference to the 
center of the cylinder. 


Now formula (42) is applied to the computation of the pressure of 
the stream on the cylinder. 


Denoting by 8 the mass density of the fluid and by X and Y 
the components of the resultant force yields, by the formula of 
Chaplygin, 


pte (ed 
c 
aW T 1 1 Ig c e 
dz oč" Zabel N i) aaa? Ne zor 


Substituting the preceding in the formula of Chaplygin and applying the 
theory of residues yield 





2gA 
X = ger? ace 
N 
diši k- kui (43) 
gr 
2 ga ži oe z cos = + A sin 5 
ST g a Go t+ de 
ne? O 288 0 
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The first formula determines the wave resistance of the cylinder 
and the second, the lift force. 


For A= - e the second formula gives the theorem of Joukowski. 
By formula (4), section 2, it is now possible to find the form of 


the free surface of the fluid. The equation of the surface of the fluid 
may be written in terms of the function W in the following form: 


c dW 
z — Reel( — + c 
mes E Ma 


Making use of formula (42), 
X z cos < (x-z) - A sin E (x-z) 
c cĉ 


T 
yolk aš 
KC n BASG 
or 
Ga E cos = + A sin či 
xC gx c c 
T y = cos > = dš + 
c -0 ECHN 
č sin g5 - A cos = 
sig | zoo. o Sa 
Co. J -o ETA 
Making x approach » gives 
EA 
S se sin & 
Y= c ae 


Hence, far behind the cylinder the surface of the fluid carries steady 
waves of length L = 2xc@/g. 


The amplitude of these waves is 


> 


g 


A = X es 
c 





f 
x 
| 
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For determining the wave resistance, energy considerations are applied. 
The wave resistance R is determined in terms of the energy E of the 
steady wave by the formula 


where U is the group velocity, and 


The energy E is determined by the formula 
E = 5 ea" 


Making use of all these results formula (43) is again found: 


9. MOTION OF CYLINDER UNDER SURFACE OF 
FLUID OF FINITE DEPTH 


Under consideration is a circular cylinder of radius R moving 
with constant velocity c under the surface of a fluid of depth H. 
The problem is proposed of finding its wave resistance as a function 
of the submersion depth h of the center of the cylinder. This prob- 
lem is solved under the assumptions of Lamb; that is, the magnitude 
h/R is assumed to be small and the entire cylinder is replaced by a 


dipole of moment -2ncRŽ, 


The motion of the fluid is studied in relation to a system of 
axes of coordinates x and y  displaced uniformly together with the 
cylinder. The velocity of the approaching flow at infinity is c. 

By (x,y) is denoted the potential of the wave velocities. The har- 
monic function %(x,y) satisfies for y= - H the condition 


O? _ 
577" (1) 
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Along the free boundary y = O the function (x,y) satisfies the fol- 
lowing condition: 


32 d$ g O$ 
NE: Kp oF O (2) 


The expression p >O is the Rayleigh coefficient of the dissipating 
forces. 


Near the point y= - h, that is, near the center of the cylinder, 
the function (x,y) on the basis of the above considerations must have 


the following form: 


cRÊx 


- ———t. 
x2+(y+h)?2 
Then 


N cRŽx (3) 


and 
2 


R x 
o sa S Rx (4) 
E x@+(y-h)? 


where % is the velocity potential of a certain fictitious dipole 
situated at the point y= h. 


In place of the function © a new function (x,y) is introduced, 
which is valid for the entire region occupied by the fluid, setting 


$ = $ -t + O (5) 


The harmonic function (x,y) satisfies the following conditions on 
the boundary of the fluid: 


a sa EN ae ' 

dy dy dy for y H (1*) 
do oo , g do | | 2g dta W : 
a ths + W; a W F for y=0 (2') 
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The function (x,y) is obtained in the form of the following definite 
integral: 


(A cos kx+C sin kx) cosh ky dk + | (B cos kx+D sin kx) sinh ky dk 


O O 

(6) 
where A, B, C, and D are unknown functions of the variable param- 
eter k. To determine these functions condition (2') is first considered: 


oF 2 E 
= E = 2 cosh Rx i = cRŽ ke = sin kx dk 
BA y=0 (x +h ) O 


Substituting this result and (6) in condition (2') yields 


œ 


( kA + uke + 85 2)008 kx + kC - Uka + S 0) sin ax 
o c c 


= - £8 p? ke ak sin kx dk 
Z o 


This relation gives the following two eguations: 


2 
c 
2 2 hk (1) 
S — 48 s 
HA + 10 - 5 D= ŠA 
Employing condition (1') gives 
%1 2cRĉx(H-h) 2 -(H-h)k 
S = - E2, vaja = CR ke sin kx dk 
Y Jyoop [x +(H-h) ] O 
2 

o? = Zon zlih) = + cR? ke (H+h)k sin kx dk 
oy y=-H [x +(H+h) ] O 
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Applying the results to condition (1') yields 


f fez cosh kH - A sinh kE) cos kx+k(D cosh kH -C sinh kH) sin 1x| a 
0 


00 


= 2cR° ke"fK sinh kh sin kx dk 


from which 
B cosh kH - A sinh kH = O 


D cosh kH - C sinh kH = 2cR e W sinh kh 


or 


B = A tanh kH 
(8) 
-Hk sinh kh 


= 2 
D = C tanh KH + 2cR" e cosh kE 


These equations together with equations (7) permit determining the 
unknown functions A, B, C, and D. 


From equations (7) and (8) the two equations for A and C are 
found: 


Ev: 


( - & tanh sr) A - jC = 0 
Cc 


g _ 28 „2 cosh k(H-h) 
HA se ie tanh sejo a Me uv 


Solving these equations determines A and C 
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28 RŽ cosh k(H-h) 


c cosh kH 
MO Fetes Stem me 
ue + (k - £ tanh kH) 


2g p2 cosh k(H-h) 
cosh kH 


u2 + (k - S tanh kE)? 
C 


(x - 5 tanh ki) 


Equations (8) then permit finding B and D. In this way the func- 
tion © is determined, which makes it possible to obtain the velocity 
potential &. The complex potential of the absolute velocities must 
be found: 


w= 6+ 10 


First 








RŽ 
o = W 2) 
p = Reel [(an1) cos kz + (C-Bi) sin = dk 
O 

where z = x + iy. Hence 

2cRĉhi 

Waon + (A+D1) cos kz + (C-Bi) sin kz] ax 
2°+h O 


The complex potential W of the relative'velocities is obtained by 
adding to the right side the term -cz. In terms of this potential 
the pressure of the stream on the cylinder is obtained by the formula 


of Chaplygin: 
2 
sis CU 
x - av = be) (22) az 


where the integral is taken over a contour containing the dipole 
z= - ih. 





30 NACA TM 1335 


10. COMPUTATION OF WAVE RESISTANCE 


The problem of finding wave resistance is restricted only to the 
computation of X, and therefore only to the imaginary part of the inte- 
gral of the formula of Chaplygin: 


2 
1 daw 
X = - Z 8 Imag. (GE) dz (9) 


The function W may be written 


cr? 
W= =- a + O2) 
The function G(z) is holomorphic about the point z = - hi: 





2 fo a] 
G(z) = - cz + PSN JI [(a+01) cos kz + (C-Bi) sin oč (10) 
O 


From formula (9) is obtained 


from which 


2 2 
(=) az = cn G'(z) dz = 4xicR"G"(-ih) 
dz (zrni) 


Therefore 
X = - 2nR°Sc + Reel G"(-ih) 


Making use of formula (10) the real part of the second derivative of 
the function G(z) at z < - ih is determined: 


Reet [o"(-sn)] = na k*(A cosh kh - B sinh kh)dk 
O 





= 
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from which 


X = 21R?8c «(A cosh 


O 


But, as was found in formula (8), B = A tanh kh. 


vo 


X = 2nRSc kĉA 


O 


Replacing A by its value gives 


cosh? k(H-h) 


X = 416gR*u 5 
cosh” kH 


O 
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kh - B sinh kh) dk 


Hence 


cosh k(H-h) dk 


cosh kH 
Kak 

u + (k - & tenh kH)? 
Cc 


which is the formula for the wave resistance in the presence of dis- 


sipative Rayleigh forces. 


cosh? k(H-h) | 
cosh® kH 


L = lim 


j=0 


H . 
O 


A new variable of integration č is 
with the following notations: 


gH 
= = x 
ce 
a 
H 
HE = Hj 


In this notation L is rewritten 





Freed from these forces 


kĉdk 


u? + (k - 5 tanh km)? 
Cc 


introduced, setting kH = č and 
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2 2 
L= lim E. fn) 
H° ugso Jo cosh” & Hy *+(& - x tanh 8) 


If the number x is less than unity 


c > gi 


then L =O and the wave resistance is equal to zero. When x> 1 
the roots of the equation in E are considered: 


2 
uy” + (E - x tanh E)? = 0 (12) 
This equation may be resolved into the two following equations: 


E 


il 


x tanh š tuji 


E 


x tanh € - uji 


Only the first of these equations will be considered. In view of the 
fact that x > 1, the equation = x tanh šč has one real root which 
is denoted by čo. Equation (11) has one real root &%, approaching čo 


as Hj approaches zero; & = a E. By applying these notations, the 
H= 


left side of equation (12) is expanded in a series about the point E: 


na? + (€ - x tanh k)? |e -x tanh & - ui | -x tanh E+ uni 


d 
= 2uji JE E -x tanh š-uji < (g-E) +... (13) 
ši E=E 


= 2u,ilf1 - t) JE . 
nif cosh E (€ z) + 


The integral (11) may be represented as the sum of a certain contour 
integral taken along the path OlK» and a residue multiplied by 2ri 
of the function under the integral relative to the point & = &: 


ta p cosh%(1-a)€ ; E | 


L= L lim + lim a 5 Z 5 
cosh" & uj(g - x tanh €) 


R 
p? pj=0 J(OIK) — uj<o E" &=$ 
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The first component on the right side gives zero in the limit. From 
series (13) the limit of the second component can likewise easily be 
found. The following expression is finally found: 

E x čo? cosh%(1-a)&, 
Hê cosh? Eo -X 


For the expression for X, 


2 
x = 425 gRX . čo cosh? (1-0) Ep 


(14) 

n? coshĉčo -x 
is now obtained. This is the formula for the computation of the wave 
resistance for velocities c less than the critical velocity vei. 
The value & is the real root of the equation 


čo = X tanh čo (15) 


From formula (14) is readily obtained the formula of Lamb for the 
wave resistance of a cylinder moving under the surface of a fluid of 
infinite depth. When H is very large, the root Eo of equation (15) 
is čo =X, from which 


zo? cosh? (1-a)zo 2.20% 
M ——————— en e 
cosh? čo-x 


Substituting this expression in formula (14), the formula of Lamb is 
obtained: 


2gh 
4 ` 2 


x, = ane (2) e f 


Now formula (14) is investigated. In place of čo is written šč. In 
the second factor on the right side x is replaced by č/tanh —E to 
obtain 


_ 8n?6gR* , £? tanh č : cosh? (1-a)č 


ae sinh ZE - ZE (18) 


X 
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The parameter x varied from 1 to œ so that & will vary from zero 
to infinity. For &=0 the second factor on the right side of this 
formula has a value equal to 3/4. It will be shown that this value is 
a minimum or maximum. Denoting by A the investigated factor gives 


Mea, z sinh 2z | [gost 0008] 
cosh & 


Then 


2 


from which 


a= $E +% - ate-a) Ea N +] 


which shows that 


ARN 
me 
knn 
o 
u 
o 


Thus for &= O the wave resistance will be an extreme value. It will 
be a minimum if a< L - 8/15; for larger values of a the vave 
resistance is a maximum. 


For E = the factor A becomes zero. This investigation 
explains the general features of the dependence of the wave resistance 
on the velocity. 
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The parameter x varied from l to e so that & will vary from zero 
to infinity. For E=0O the second factor on the right side of this 

formula has a value egual to 3/4. It will be shown that this value is 
a minimum or maximum. Denoting by A the investigated factor gives 


poi Ei sinh 2E [sään 008] 
-i, . [cosh (1-a)č 


sinh 2k - 2č cosh č 
Then 
E E er Ma PR 
Z "sim ŽE ŽE 4 15 oy Pe va 


2 
[esse Gnade)” = a- aea? + 


from which 


a~ ip + (1 - at-a) ore id 


which shows that 


ARN 
R 
uN, 
o 
li 
o 


Thus for &= O the wave resistance will be an extreme value. It will 
be a minimum if a< 1 - 48/15; for larger values of a the wave 
resistance is a maximum. 


For € =” the factor A becomes zero. This investigation 
explains the general features of the dependence of the wave resistance 
on the velocity. 
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11. CONCLUSIONS 
Formula (16) obtained in section 10 for the wave resistance of a 
cylinder of radius R moving under the surface of a fluid of finite 
depth H led to the following conclusions with regard to the change in 
the wave resistance with velocity of motion and with depth h of its 
submersion: 


The magnitude 


was computed for a series of values of a= h/H: 


3 1 8 0327 15 
3 Ig’ 7? 1- affs 0.3; 0.36; 0.4; 0.5; g? 7’ a? 16 


for the Froude number f = c/ yg varying between zero and 1. The 
results of the numerical computations are given in figure 1. 


e 
li 
œj = 


The curves of this diagram show a considerable wave resistance 
for small depths of submersion of the cylinder. With increased sub- 
mersion of the cylinder the wave resistance drops sharply for most of 
the values of the Froude number. 


For the value of the parameter a=1 - 8/15, the character of 
the maximum or minimum changes at the point f= 1. For values of 
a< 1 - 8/15 the curves of the wave resistance have a sharply formeä 
maximum. With a increasing from zero to 1 - ,/8/15 this maximum 
decreases. The presence of this maximum shows up also on the curves 
corresponding to values of a somewhat less than the number 1 - 4/8/15; 
the curves of wave resistance for a between 1 - V8/15 and 0.3 
(approximately) have two peaks, one at the point f < l and the other 
near the point f = 0.55. For values of a> 0.3 the wave resistance 
increases monotonically with increasing velocity; the rate of increase 
of the wave resistance is considerable for Froude numbers f near the 
critical number f =1. 


In regard to the problem of the first part of the paper (sections 2 
to 8), the results may be described as follows: When the motion of a 
circular cylinder under the surface of a fluid is studied, the cylinder 
may be replaced by a vortex if the ratio of the radius of the cylinder 
to the depth of submersion of its center is less than 1/3 (see table). 
In the composition of the flow and the lift force the principal part is 
played by the circulation of the velocity about the cylinder. However, 
by considering the general solution of the problem the terms obtained 





A ==. 
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by Lamb and Havelock in their work are found. The effect of these 
terms is appreciable only for small depths of submersion of the 
cylinder; for larger depths the circulation of the stream velocity 
plays the fundamental part. 


Translated by S. Reiss 
National Aävisory Committee 
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